Binary-Induced Gravitational Collapse: A Trivial Example 
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We present a simple model illustrating how a highly relativistic, compact object which is stable in 
isolation can be driven dynamically unstable by the tidal field of a binary companion. Our compact 
object consists of a test-particle in a relativistic orbit about a black hole; the binary companion is 
a distant point mass. Our example is presented in light of mounting theoretical opposition to the 
possibility that sufficiently massive, binary neutron stars inspiraling from large distance can collapse 
to form black holes prior to merger. Our strong-field model suggests that first order post-Newtonian 
treatments of binaries, and stability analyses of binary equilibria based on orbit-averaged, mean 
gravitational fields, may not be adequate to rule out this possibility. 



I. INTRODUCTION 

Binary neutron stars are among the most promising 
sources for gravitational wave detectors now under con- 
struction, like LIGO, VIRGO and GEO. This fact has 
motivated an intense theoretical effort to understand 
their dynamical evolution and predict the gravitational 
waveforms emitted during their inspiral and coalescence. 

Fully general relativistic treatments of the binary prob- 
lem are nontrivial because of the nonlinearity of Ein- 
stein's field equations and the need for very large com- 
putational resources to solve them, together with the 
equations of relativistic hydrodynamics, in three spatial 
dimensions plus time. Recently, Wilson, Mathews and 
Maronetti (hereafter WWM; performed approximate 
relativistic simulations which suggest that neutron stars 
just below the maximum allowed rest mass and dynam- 
ically stable in isolation become unstable to collapse to 
black holes when placed in a close binary orbit. Their re- 
sults are in disagreement with earlier calculations of sta- 
bility performed in Newtonian gravitation, which show 
that tidal fields stabilize binary stars against radial col- 
lapse Q . The WWM results are also surprising in light 
of several recent, albeit simplified, post-Newtonian (PN) 
perturbation analyses || jl| , none of which indicate the 
presence of any relativistic radial instability in binaries, 
at least to first post-Newtonian (1PN) order. 

To help clarify the issue, Baumgarte et. al. || per- 
formed fully relativistic numerical calculations of binary 
neutron stars in quasi-equilibrium circular orbits. (The 
solutions correspond to "quasi" -equilibrium because of 
the inevitable slow, secular decay of the orbit due to the 
emission of gravitational waves). They found that the 
maximum allowed rest mass (i.e. baryon mass) of a neu- 
tron star in a synchronous binary is in fact slightly larger 



than in isolation. PN treatments of the same problem 
confirm this conclusion and extend it to nonsynchronous 
binaries as well (e.g. Ref [||). This finding rules out at 
least one explanation for the results of WWM - namely, 
the possibility that binary equilibrium solutions might 
not exist for neutron stars with the highest values of rest 
mass allowed for isolated stars. Using their relativistic 
models, Baumgarte et. al. [|| further showed that all 
stars which are stable to radial collapse in isolation are 
also stable when placed in a binary, at least to secular 
modes, all the way down to the innermost stable circu- 
lar orbit. This result, together with the expectation that 
dynamical instabilities in binaries always arise at smaller 
separation than secular instabilities, as in Newtonian the- 
ory @] ' suggests that a dynamical instability against 
radial collapse is just not present for binary neutron stars 
that inspiral from large distance. This conclusion has 
been considerably strengthened by an analytic, relativis- 
tic analysis by Thorne g] , who argues that a rigidly ro- 
tating, stable neutron star in an inspiraling binary system 
will be protected against secular and dynamical collapse 
by tidal interaction with its companion. 

It would appear that the WWM result may not be 
correct and that their finding might be the result of in- 
adequate computational resolution or one or more of the 
simplifications they adopt to solve this difficult numerical 
problem. Nevertheless, their work has raised the very in- 
teresting and nonobvious question about whether or not 
tidal fields in a relativistic binary can under any circum- 
stances trigger the collapse of a compact object that is 
known to be stable in isolation. 

In this paper, we offer a simple example to illustrate 
the possibility of binary-induced dynamical collapse of a 
compact object. The compact object consists of a test 
particle in a tight orbit around a black hole. This "com- 
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pact object" is itself placed in orbit about a distant mass 
and the "binary" - i.e., the "compact object" in circu- 
lar orbit about the distant mass - is allowed to inspiral 
by the emission of gravitational radiation. We investi- 
gate whether the tidal field of the distant mass can drive 
such a compact object dynamically unstable even though 
the object, when in isolation, is stable. Here, dynamical 
instability manifests itself as a rapid plunge of the test 
particle toward the black hole. 

Our analysis of this model is highly simplified and ar- 
guably heuristic in parts. Nevertheless, we are confident 
that our overall picture is reliable and that our essential 
conclusions will stand up to more rigorous scrutiny. 



II. BASIC MODEL AND EQUATIONS 

Consider at first a self-gravitating, N-body system 
in the post-Newtonian limit of general relativity. The 
(geodesic) equations of motion for body i may be cast in 
the form 
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where rrij is the total mass-energy of body j, rj is its 
coordinate position and ry = rj — rj . The first term on 
the right hand side of Eq. ([!]) is the usual Newtonian 
(Coulomb) expression, while the second term is the PN 
correction. Explicit expressions for this PN correction 
have been developed by many authors to various orders 
and in several gauges. Einstein, Infield and Hoffman ||] 
provided the "EIH" equations of motion, which are valid 
to 1PN order. Damour and Shafer pi derived an N- 
particle Lagrangian valid to 2PN order (see Ref for a 
review and references). In these treatments, the PN cor- 
rection to the equation of motion for a particle typically 
depends on the positions, velocities and the accelerations 
of all of the N particles. 

Let us now specialize to TV — 3. Assume that bodies 
1 and 2 are appreciably closer to each other than they 
are to body 3, and treat the influence of 3 on the relative 
orbit of the close pair 1 — 2 as a small perturbation. Define 
the quantities m = mi + m-i and 

miri + m 2 r 2 _ 
Rem = , r = n-r 2 , r cm = R cm -r 3 , 

m 

(2) 

and use Eq. (^) to write the equation of motion for the 
relative orbit of the close pair in the form 



is the leading Newtonian tidal perturbation due to body 
3 and where the leading post-Newtonian corrections Fi 
and F 2 depend only on bodies 1 and 2. Here we assume 
that r <C r cm and we drop tidal corrections of 0{m^^^). 
To the same order, the orbit of particle 3 about the close 
pair is determined by 



M 



(5) 



where M = 
G(M\ 



m + m 3 and where we drop tidal terms of 
as well as the very small PN corrections. 

Here m> = max[mi, 7712], and similarly for m < . To this 
order, the orbit of body 3 about the close pair decouples 
from the relative orbit of that pair. 

In this perturbative treatment, let us henceforth ex- 
amine two extreme, opposite limits: In one limit, assume 
that the orbital plane of body 3 about the close pair 
moves in the plane describing the relative orbit of that 
pair (i.e., the inclination angle of the two orbital planes 
satisfies i = 0°). In the other limit, assume that 3 moves 
in a plane perpendicular to the orbital plane of the pair 
(i = 90°). Decompose Eqs. (||) and (||) into radial and 
angular components, yielding for the relative orbit of the 
pair 
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and for the orbit of 3 about the pair 
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In the Newtonian limit, A — 1 and B = in Eqs. (g) 
and ( f7[). In the case of isolated binaries, Lincoln and 
Will ]12|] derive post-Newtonian expressions for A and 
B for arbitrary masses, correct through 2.5PN order. 
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Kidder, Will and Wiseman (hereafter KWW; |13|]) pro- 
vide a "hybrid" set of equations in which the sum of 
the terms in A and B that are independent of the ratio 
r\ = (Jt/m, fi = m\mi/m^ is replaced by the exact expres- 
sion for geodesic motion in the Schwarzschild metric of a 
body of mass m, while the terms dependent on rj are left 
unaffected. Their resulting equation of motion is there- 
fore exact in the test-body limit (rj — > 0) and is valid to 
2.5PN order when appropriately expanded for arbitrary 
masses. We shall utilize these same hybrid expressions 
for A and B and, for simplicity, work in the test-body 
limit by taking one member of our close pair to have a 
mass much smaller than the other. Adopting harmonic 
(or de Donder) coordinates, the resulting (Schwarzschild) 
expressions for A and B are given by 
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We are interested in the dynamical behavior of the 
close pair, regarded as a single "compact object", as it 
inspirals toward the distant mass 7713 due to gravitational 
radiation emission. To treat the inspiral of this "binary" 
(7773 in orbit about the "compact object"), we must in- 
clude radiation reaction terms in the lowest order (New- 
tonian) orbit Eqs. ( |l0| ) and (fTl]). Formally, such a treat- 
ment requires a consistent expansion up to 2.5PN order. 
In lieu of this, we shall analyze the inspiral by assuming 
that the binary is in a nearly circular, Keplerian orbit, 
which undergoes a slow inspiral due to gravitational ra- 
diation loss in the quadrupole limit. This assumption 
is equivalent to inserting a quadrupole radiation reac- 
tion potential in the binary orbit Eq. (g) and neglecting 
the lower-order, (non-dissipative) PN corrections in that 
particular equation. While such an expression is not for- 
mally consistent to 2.5PN order, we believe that it faith- 
fully tracts the secular inspiral of the binary in the limit 
treated here in which the binary system is at wide (non- 
relativistic) separation. The details of the inspiral are 
not important here, only that the inspiral serves to bring 
a tidal perturber slowly in from infinity toward our com- 
pact object [[l4). The resulting equation for the binary 
inspiral is then 



r cm (t)/r cm (0) = (1 - t/T) 1 / 4 
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where the binary inspiral timescale T is given by 
5 (r cm (0)/m) 4 
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We will use Eqs. (|1 5|) — (pL7[) in analyzing the relative orbit 
Eqs. (§)-(!). 

All orbital timescales in our problem are considerably 
shorter than T. For some purposes it is instructive to 
perform an orbit average of the tidal terms appearing in 
Eq. (g) and (0), which yields {f( idal ) = am 3 r/r^ m , where 

a = 1/2 for i = 0°, -1/4 for % = 90°, and (/ t f dal ) = 0. 



III. QUASI-EQUILIBRIUM AND STABILITY: 
MEAN-FIELD APPROACH 

Use the orbit-averaged tidal force to consider the close 
pair as a compact object in the mean tidal field of its dis- 
tant, binary companion. In this approximation the mean 
field, quasi-equilibrium state for the compact object ob- 
ject consists of a test particle in circular orbit at constant 
radius about a black hole. The degree of compaction of 
this quasi-equilibrium compact object is parametrized by 
the radius of the circular orbit, ro/m; for small values of 
this parameter < 10, the object is relativistic. The quasi- 
equilibrium orbit in the mean tidal field is determined 
self-consistently by requiring r = = f = (j>. Eq. (|^) 
then implies the circular orbit condition 
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which can be solved along with Eqs. (|T^) and (|l4|) to get 
f^o as a function of rQ. Because the orbiting particle is a 
test-particle of negligible mass, its orbit does not decay 
due to gravitational radiation, so the compact object is 
secularly stable. Dynamical stability of the circular orbit 
is determined by considering linear perturbations about 
the equilibrium values of ro, 4>o and r'o = 0. The result- 
ing stability analysis of our compact object is a simple 
extension of the analysis by KWW for PN binaries in the 
absence of a tidal field. We find by linearizing Eq. (g) in 
a mean tidal field that the condition for stable circular 
orbits may again be written as a + be < 0, where now 
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We can analyze some limiting cases analytically. In the 
absence of a tidal field (r cm /m — > 00), the condition for 
dynamical stability reduces to 
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FIG. 1. The innermost stable circular orbit (ISCO) as a 
function of the separation between the compact object and its 
binary companion in the mean-field approximation. The solid 
line is for coplanar test-particle and binary companion orbits 
(i — 0°); the dashed line is for orthogonal orbits (i — 90°). 
The dotted line shows the ISCO for an isolated compact ob- 
ject. 

Eq. (^o|) reproduces the well-known value (in harmonic 
coodinates) of the innermost stable circular orbit (ISCO) 
in Schwarzschild geometry, r$/m — 5. Stability for a 
strictly Newtonian orbit (ro/m ^> 1) in the presence of a 
tidal field from a perturber in the orbital plane requires 

r /r cm < (- -) , fig > ^ (21) 

(Newtonian, i = 0°), 

which is the familiar Newtonian condition for avoiding 
tidal disruption. For a tidal perturber perpendicular to 
the orbital plane, all Newtonian orbits are stable. 

The more interesting cases involve the effect of the 
mean tidal field of a companion on the stability of a rela- 
tivistic compact object in binary equilibrium. To analyze 
these, we numerically evaluate Eq. ( |l9| ) to determine the 
ISCO as a function of the strength of the tidal field; re- 
sults are summarized in Fig 1. We set m — so that 
the compact object and its distant companion have equal 
mass. We see that for a perturber in the orbital plane of 
the close pair, the ISCO moves out as the perturber ap- 
proaches and the tidal field increases. Sufficiently tight 
orbits with r^jm > 5, which are always stable when the 
compact object in isolation, are destabilized by the com- 
panion. To guarantee that this result is not just a co- 
ordinate effect, we can parametrize an orbit not only by 
its radius, which is not gauge- invariant, but also by its 
period or angular frequency Sl n , which can be measured 



FIG. 2. Angular frequency at the ISCO as a function 
of the separation between the compact object and its binary 
companion in the mean-field approximation. Curves are la- 
beled as in Fig 1. 

by a distant observer. We see from Fig 2 that for a per- 
turber in the orbital plane, the angular frequency of the 
ISCO moves to lower frequencies with increasing tidal 
field strength. Thus, a given compact object, specified 
by a fixed angular frequency, can be stable in isolation 
but driven dynamically unstable when in a binary sys- 
tem. 

From Fig 1 we see, by contrast, that for a perturber 
perpendicular to the orbital plane the ISCO actually 
moves in. Likewise, we see from Fig 2 that the orbital 
frequency at the ISCO increases. Hence, in this case the 
mean tidal field due to a companion would appear to 
stabilize the compact object. 

IV. STABILITY DURING INSPIRAL 

To assess the true fate of our compact object inspiral- 
ing toward its binary companion, we abandon the mean- 
field approximation and integrate the relative orbit Eqs. 
(|)-(|) together with the inspiral Eqs. (|l|)- (|l7|). At 
t = we place a test-particle in a circular orbit about 
a black hole, setting r/m = 5.9, r = = <p an d </> given 
by Eq. (|l8|). We set 7713 = m and start the companion 
ms at large separation r cm /m = 40. Hence, initially, the 
tidal field of the companion is negligible at the compact 
object and with so large a test particle orbit, the com- 
pact object begins in stable equilibrium, almost exactly 
as it would be were it in isolation. We again treat two 
cases for the relative orientation of the two orbital planes, 
i = 0° and i = 90°, setting 9 = 0° at t = in the former 
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FIG. 3. The orbital radius of the test particle about the 
black hole as a function of time. The solid curve shows the 
motion for coplanar test-particle and binary companion orbits 
(i = 0°); the dashed curve shows the motion for orthogonal 
orbits (i — 90°). The dotted curve shows the separation be- 
tween the compact object and its binary companion as they 
inspiral together. 

case (companion aligned with the close pair) and 9 = 90° 
in the later case, (companion situated above the orbital 
plane of the close pair, along its normal) At such large bi- 
nary separation, QqT 1, hence the binary orbits many 
times before the inspiral proceeds very far. We thus ex- 
pect (and have verified numerically) that the fate of the 
compact object does not depend on our choices of initial 
values of 4> and 6. 

The results of our integrations are summarized in Figs 
3-5, where we see that early on the compact object re- 
mains in stable equilibrium for many binary periods as 
the companion slowly approaches. The test-particle or- 
bit during this phase is characterized by stable, small 
oscillations about the initial circular orbit; the oscilla- 
tions occur at the beat frequency between the distant 
binary companion and test-particle orbits (cf. Eq (pD). 
When the companion gets sufficiently close (r cm /m ss 20 
for i = 0°, r cm /m w 30 for i = 90°), the amplitude of 
the oscillations becomes sufficiently large to cause the 
test-particle to plunge toward the hole and be captured. 
Hence in both cases the tidal field of the companion ul- 
timately drives the dynamical collapse of the compact 
object. 

When i = 0°, this binary- induced collapse is antici- 
pated from the mean-field stability analysis, since in this 
case, the tidal field destabilizes close circular test-particle 
orbits which are otherwise be stable in isolation (see Sec 
3). The unexpected destabilization in the case i = 90° 
appears to result result from the growing eccentricity 
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FIG. 4. The trajectory of the test particle about the black 
hole during binary inspiral for a companion moving in a copla- 
nar orbit (i = 0°). 
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FIG. 5. The trajectory of the test particle about the black 
hole during binary inspiral for a companion moving in an 
orthogonal plane (i — 90°). 
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of the test-particle orbit as the companion approaches, 
eventually driving the particle too close to the hole at 
pericenter and to ultimate capture. In this case, the 
mean- field analysis, based on orbit averaging and strict 
circular motion, does not account for the more compli- 
cated motion of the test particle and fails to predict its 
ultimate fate. 



V. DISCUSSION 

We have presented a simple illustration of a compact 
object which is stable in isolation but dynamically un- 
stable to collapse when inserted in a binary. Admittedly, 
our compact object is highly idealized and our treatment 
is simplified. Our whole system consist of only 3 bodies 
moving on geodesies. We could trivially generalize our 
compact object to be a fully three-dimensional swarm 
of test particles, all in randomly oriented circular orbits 
at the same fixed radius about a central hole. When- 
ever this radius is close to, but slightly larger than 5M, 
then according to our previous analysis, the presence of 
a sufficiently close binary companion would drive most, 
if not all, of the test-particle orbits dynamically unstable 
during inspiral from large separation. Once again, such 
a compact object would be stable in isolation but could 
be driven dynamically unstable when placed in a binary 
system. 

It is not clear what bearing our simple example has to 
the finding of WWM, especially in light of the counter ev- 
idence suggesting that binary-induced collapse does not 
occur in fluid stars. We have not incorported any hy- 
drodynamic forces in our model, and these forces are 
necessary to prevent catastrophic collapse in fluid sys- 
tems. It is interesting nevertheless that a simple colli- 
sionless particle model exists which is stable in isolation 
but undergoes dynamical collapse during binary inspiral. 
Typically, collisionless equilibria in isolation are subject 
to the many of the same relativistic instabilities as fluid 
systems, including radial instability to collapse, (see, eg, 
Ref |15| for a review and references), so it is certainly 
relevant for WWM and others to have considered the pos- 
sibility of binary-induced collapse in the case of binary 
neutron stars. 

Our simple model calculation offers a few cautions re- 
garding the treatment of the relativistic binary problem. 
Crucial to our finding was our use of a high order rel- 
ativistic equation for the dynamical behavior of of our 
compact object. Indeed, had we used 1PN expressions 
for A and B in Eq s. (o ) and (0) instead of the hybrid 
expressions Eqs. (OF and (OJ we would have found 
spurious behavior, because at 1PN (and even 3PN), the 
solution for the test-particle ISCO in a PN expansion of 
the Schwarzshild equations of motion yields a spurious 
root that does not converge to 5m for isolated holes [OJ . 
By using the hybrid expressions in Eqs. @ and (0), our 
compact object is treated exactly to all orders, at least 



when it is isolation. 

It is not surprising that a 1PN treatment is inadequate 
for this problem. Such an approximation is not even suf- 
ficient to determine reliably the onset of radial instability 
to collapse in an isolated fluid star. There the highly non- 
linear terms of Einstein's equations play a crucial role and 
only for special equations of state (e.g n = 3 polytropes) 
is a 1PN analysis quantitatively reliable Q|. 

We have also seen that stability analyses of equilibria 
based on orbit averaging and mean tidal fields may prove 
misleading, as they did for the i = 90° case treated above. 
By the same token, it is by no means clear that conclu- 
sions found for synchronous stellar binaries will apply to 
nonsynchronous binaries. (Here we recall that the vis- 
cosity in neutron star binaries is not sufficiently large to 
drive the system to corotation jl6]]). As our calculations 
demonstrate, once we relax dynamical constraints and 
allow motion with greater degrees of freedom, the final 
outcome of binary inspiral can be quite different, even 
qualitatively. 

There is mounting evidence which argues against the 
possibility that sufficiently massive, highly compact neu- 
tron stars in coalescing binaries can collapse to black 
holes prior to merger. However, we will be greatly re- 
assured when this conclusion is finally corroborated by 
detailed hydrodynamical calculations in full general rel- 
ativity. 
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